one_way_res = [0 for _ in range(16)]
rnd_perm = [[0 for _ in range(16)] for _ in range(256)]
rnd_glob = [0 for _ in range(256)]
ciphertext = [0 for _ in range(32)]
compt, i, y = 0, 0, 0
l, r = [[[0 for _ in range(16)] for _ in range(32)]] * 2

chaos = [0 for _ in range(256)]

chaos1 = [173,84,240,67,1,53,254,36,41,172,115,109,223,199,152,189,90,46,
149,193,218,130,250,40,203,4,35,237,236,246,213,143,169,176,48,
23,61,206,69,34,97,155,4,109,183,220,42,64,21,123,29,233,253,
105,183,209,1,191,113,12,46,7,8,183,166,199,166,7,78,37,135,
252,174,84,140,164,152,94,22,185,59,68,181,60,176,67,51,25,28,
190,138,198,44,90,92,221,149,175,186,25,49,210,50,237,41,207,
31,226,114,121,230,15,58,25,142,58,98,232,59,3,189,28,8,116,
131,185,78,250,239,33,116,173,94,45,104,62,122,179,18,150,246,
250,17,8,79,157,225,238,47,10,133,58,8,126,82,68,153,141,2,158,
204,50,130,53,59,32,243,160,172,35,24,107,35,115,228,143,28,
224,77,55,25,28,120,89,186,152,49,84,117,180,30,138,134,77,182,
157,61,230,22,149,54,15,110,32,213,155,106,78,16,23,89,140,158,
169,96,136,186,104,30,199,67,35,218,159,210,109,28,238,33,150,
173,180,247,201,83,150,105,164,228,59,207,101,221,99,52,120,
199,31,6,144,202,215,209,49,42,195]


# Not Done
def init_pcb1(key) -> int:
    global one_way_res, rnd_perm, rnd_glob, ciphertext, compt, l, r, chaos, chaos1
    
    w, err, ix, x = 0, 0, 0, 0
    
    fixed_key = [46,245,138,13,244,233,238,154,139,30,241,90,47,205,171,97,223,
  190,28,10,185,13,23,137,30,208,254,143,165,101,27,48]

    for compt in range(0, 256):
        chaos[compt] = chaos1[compt] ^ key[compt%20]
    
    for w in range(0, 4):
        round_init(fixed_key, 32)

        for compt in range(0, 256):
            chaos[compt] = chaos[compt] ^ ciphertext[compt%32]

        for compt in range(0, 32):
            fixed_key[compt] = ciphertext[compt]
    
    for compt in range(0, 16):
        err = 0 
        for ix in range(0, 16):
            round_init(ciphertext, 4)

            for x in range(0, ix):
                if rnd_perm[compt][x] == (ciphertext[7]%16):
                    x = ix
                    err = 1
                    round_init(ciphertext, 4)

            if err == 0:
                rnd_perm[compt][ix] = ciphertext[7]  % 16
            else:
                err = 0
                ix = ix - 1
    
    err = 0
    for ix in range(0, 256):
        round_init(ciphertext, 4)

        for x in range(0, ix):
            if rnd_glob[x] == (ciphertext[7] % 256):
                x = ix 
                err = 1
                round_init(ciphertext, 4)
        
        if err == 0:
            rnd_glob[ix] = ciphertext[7] % 256
        else:
            err = 0
            ix = ix - 1
    
    for compt in range(0, 256):
        chaos[compt] = chaos[compt] ^ key[compt%20]
    
    for w in range(0, 4):
        round_init(fixed_key, 32)

        for compt in range(0, 256):
            chaos[compt] = chaos[compt] ^ ciphertext[compt % 32]
        
        for compt in range(0, 32):
            fixed_key[compt] = ciphertext[compt]
            

def round_init(plaintext, nbround):
    global one_way_res, rnd_perm, rnd_glob, ciphertext, compt, l, r, chaos, chaos1

    for i in range(0, 16):
        l[0][i] = plaintext[i]
        r[0][i] = plaintext[i+16]

    for i in range(1, nbround):
        one_way_init(r[i-1])

        for y in range(0, 16):
            l[i][y] = r[i-1][y]
            r[i][y] = l[i-1][y]^one_way_res[y]
    
    for i in range(0, 16):
        ciphertext[i] = l[nbround-1][i]
        ciphertext[i+16] = r[nbround-1][i]


def round(plaintext, nbround):
    global one_way_res, rnd_perm, rnd_glob, ciphertext, compt, l, r, chaos, chaos1

    for i in range(0, 16):
        l[0][i] = plaintext[i]
        r[0][i] = plaintext[i+16]
    
    for i in range(1, nbround):
        one_way(r[i-1])

        for y in range(0, 16):
            l[i][y] = r[i-1][y]
            r[i][y] = l[i-1][y] ^ one_way_res[y]
    
    print("\n")

    for i in range(0, 16):
        ciphertext[i+16] = l[nbround-1][i]
        ciphertext[i] = r[nbround-1][i]


def one_way_init(matrix):
    global one_way_res, rnd_perm, rnd_glob, ciphertext, compt, l, r, chaos, chaos1

    z = 0
    permut = [
    [5,10,1,3,0,4,15,2,7,9,13,8,11,14,12,6],
    [5,7,15,13,3,8,1,11,10,0,9,12,2,6,14,4],
    [4,6,12,14,1,0,13,5,7,9,10,11,2,15,8,3],
    [3,4,2,12,13,9,1,8,5,15,14,6,10,0,11,7],
    [9,8,3,5,4,6,0,2,13,15,14,1,12,11,7,10],
    [14,13,3,11,0,2,1,10,7,9,4,5,15,6,12,8],
    [8,4,0,5,7,3,12,11,13,9,6,1,15,10,2,14],
    [11,14,13,1,9,8,4,5,12,6,7,3,10,0,15,2],
    [13,12,2,1,6,9,14,11,3,10,8,15,4,5,7,0],
    [1,14,12,7,8,0,11,9,6,5,4,2,10,3,15,13],
    [6,0,10,8,5,2,11,13,14,1,4,9,3,12,15,7],
    [12,0,5,13,2,4,8,6,3,15,9,10,14,7,1,11],
    [7,10,8,11,9,2,0,3,5,4,12,15,14,1,6,13],
    [8,4,15,7,10,9,3,12,6,14,0,2,1,5,13,11],
    [11,13,14,3,5,12,7,6,2,0,10,8,1,15,4,9],
    [1,6,3,8,7,4,11,14,0,10,9,2,15,12,5,13]]

    for z in range(0, 16):
        one_way_res[z] = ((chaos[matrix[permut[z][0]]] | chaos[matrix[permut[z][1]]]) + \
            (chaos[matrix[permut[z][2]]] | chaos[matrix[permut[z][3]]])) ^ \
            ((chaos[matrix[permut[z][4]]] + chaos[matrix[permut[z][5]]]) + \
            (chaos[matrix[permut[z][6]]] ^ chaos[matrix[permut[z][7]]]))

        one_way_res[z] = one_way_res[z] + ( ( (~chaos[matrix[permut[z][8]]] ^ chaos[matrix[permut[z][9]]]) + \
                (chaos[matrix[permut[z][10]]] & ~chaos[matrix[permut[z][11]]])) ^ \
                ((chaos[matrix[permut[z][12]]] ^ ~chaos[matrix[permut[z][13]]]) + \
                (chaos[matrix[permut[z][14]]] ^ chaos[matrix[permut[z][15]]])))


def one_way(matrix):
    global one_way_res, rnd_perm, rnd_glob, ciphertext, compt, l, r, chaos, chaos1

    z = 0
    for z in range(0, 16):
        one_way_res[z]=\
            ((chaos[matrix[rnd_perm[z][0]]] | chaos[matrix[rnd_perm[z][1]]]) + \
            ( chaos[matrix[rnd_perm[z][2]]] | chaos[matrix[rnd_perm[z][3]]])) ^ \
            ((chaos[matrix[rnd_perm[z][4]]] + chaos[matrix[rnd_perm[z][5]]]) + \
            ( chaos[matrix[rnd_perm[z][6]]] ^ chaos[matrix[rnd_perm[z][7]]]))

        one_way_res[z] = one_way_res[z] + \
            (((~chaos[matrix[rnd_perm[z][8]]] ^   chaos[matrix[rnd_perm[z][9]]]) + \
            (   chaos[matrix[rnd_perm[z][10]]] & ~chaos[matrix[rnd_perm[z][11]]]) ) ^ \
            ((  chaos[matrix[rnd_perm[z][12]]] ^ ~chaos[matrix[rnd_perm[z][13]]]) + \
            (   chaos[matrix[rnd_perm[z][14]]] ^   chaos[matrix[rnd_perm[z][15]]]) ) )

        one_way_res[z] = one_way_res[z] + rnd_glob[one_way_res[z]] + rnd_glob[z]


if __name__ == '__main__':
    tt = 0
    
    key = [
        1,0,0,0,0,0,0,0,0,0,
        0,0,0,0,0,0,0,0,0,0
    ]

    text = [
        0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
        0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
    ]

    print("Plaintext : " + ' '.join(f"{i:02x}" for i in text))

    #### ENCRYPT ####
    init_pcb1(key)
    
    round(text, 6)  # 5 rounds

    print("Ciphertext : ", end='')
    for tt in range(0, 32):
        print(f"{ciphertext[tt]} ", end='')
        text[tt] = ciphertext[tt]
    
    #### DECRYPT ####
    init_pcb1(key)
    
    round(text, 6)

    print("Plaintext : " + ' '.join(f"{i:02x}" for i in ciphertext))

    print("\n\n")


"""
key :      0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Plaintext :0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Ciphertext :B 58 52 F 6A 3C 1E 63 70 46 16 73 57 BC 68 CD BF 7C E0 8F 15 54 3A E5 59 96 55 54 83 43 E 86

key :      1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Plaintext :0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Ciphertext :D9 48 EE 9C 33 94 7A A3 F F0 1E F9 CB E2 32 50 BA 43 E1 18 47 D3 24 89 4A 63 14 E5 D B6 9E 2B

key :      0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Plaintext :1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Ciphertext :4 CF B 24 E9 A8 3A AD C1 53 6A 7C 46 E1 3D CD 8C A8 9B 3D 10 B9 DC A1 D3 F 44 A6 D5 1F 93 83

key :      1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Plaintext :1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Ciphertext :4E 43 60 1B 4F 57 D2 CA DA BB 9D E0 F5 61 3 D3 55 6B C5 DD C4 8F 48 65 16 A1 3A 15 B1 B8 19 58
"""
